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Abstract. The X-ray burster and dipper 4U1915-05 
(also known as XB1916-053) was observed 19 times for 
a total exposure of roughly 140 ks between 1996 Febru- 
ary and October with the proportional counter array 
(PC A) aboard Rossi X-ray Timing Explorer (RXTE). 
We present here the detailed timing study of its X-ray 
emission. The source was observed in both a high in- 
tensity/soft spectral and a low intensity/hard spectral 
state. The 2-20 keV luminosities are about 1.4 x 10 37 and 
3.2 x 10 36 ergs s" 1 in the two regimes respectively (as- 
suming a source distance of 9.3 kpc). We confirm that 
4U1915-05 is an atoll source based on its spectral be- 
haviour and its aperiodic variability in the frequency 
range ~ 4 x 10~ 3 - 128 Hz. 4U1915-05 displays red noise 
below ~ 1 Hz, but is especially remarkable for its quasi- 
periodic variability detected for the first time at several 
timescales. We detect low frequency quasi-periodic oscil- 



crcd X-ray bursters XTE J2123-058 and XTE J 1723-376 



(Marshall and Markwardt, 1999). Following their discov- 
eries, it has been realized that they could possibly be 
used to constrain the fundamental parameters of neutron 
stars, such as their spin periods, their masses, and their 
radii. Implications of general relativity in the interpreta- 
tion of these HFQPOs have also been d iscussed by severa l 
authors like, e.g., K aaret et al. ( |l997|) , Kluzniak ( 1998 ), 
Miller et al. ( |l998b| ). 

In most systems, twin HFQPOs are observed with a 
separation in the range 220-360 Hz. Changes in the source 
luminosity cause the twin peaks to shift simultaneously 
in a way that their separation remains approximatively 
constant. For the Z source Sco X-l and the atoll sources 
4U1608-52, 4U1728-34 and 4U1735-44, however, the fre- 
quency separation significantly decreases as the HFQPOs 



free uency QPUs (HFQPUs) between ~ 200 and 1300 Hz 
Both LFQFUs and HFQPOs frequencies positively cor- 
relate with the mass accretion rate inferred from the po- 
sition of the source on the color-color diagram. In a nar- 
row range of luminosity, we also detect twin simultaneous 
HFQPOs separated by 348 ± 12 Hz. 



frequencies increase (Van der Klis et al., 1997; Mcndcz 



latij ms ii-i Qi-u>! iK-iw.'c-i ~ 5 and »u tiz as well as high ^ al 1998c | : |Mcndez and Van dcr Klis, 199 9|; |Ford et al 



1998H ; [Psaltis et al., 1998| ). In a few X-ray bursters (e.g. 



4U1702-43, KS1731-260), which are atoll sources, coher- 
ent pulsations have been seen during type I X-ray bursts 
at frequencies equal or tw ice the frequency difference be- 
tween the tw in HFQPOs flMarkwardt et al., 1999| ; |Smith 
et al., 1997] ). In 4U1636-53 and 4U1728-34, burst oscil- 



lations have been detected at a frequency significantly 



1. I ntroduction 



greate r than the twin HFQPOs separati on ( Mendez et al. 
19984 |Mendcz and Van dcr Klis, 1999| ). 



1.1. High frequency quasi-periodic oscillations in 
low-mass X-ray binaries 

One of the greatest successes of RXTE is the discov- 
ery of HFQPOs in the X-ray emission of low-mass X- 



rayi binarioo (LMXBo). RXTE made it poooiblo by the 
co mbination of the largo effective area of the PCA, an 



Specific also to the atoll sources, HFQPOs are de- 
tected when the source is in the so-called island state 
and on the lower branch of the banana state. These two 
states are easily identified in color-color diagrams and 
there is conclusive evidence that the accretio n rate in- 
creases from the island to the banana state ( Hasingei 



and Van der Klis, 1986 ). They are associated with dif- 



excellent time resolution, an extended telemetry band- 
width and optimized observing efficiency. HFQPOs rang- 
ing from ~ 300 to ~ 1200 Hz have been reported so far 
from 21 neutron star LMXBs ( |Van der Klis, 1999| ): the 6 
Z sources, 12 known atoll sources, the unknown bursting 
source in the galactic center region, the recently discov- 



ferent types of power density spectra (PDS) at low fre- 
quencies (in the remainder of this paper we define low 
frequencies as frequencies below 100 Hz). In the island 
state, the source displays high frequency noise (HFN also 
called "flat topped noise") and the PDS can be approxi- 
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mated by a zero-centered Lorentzian ( Olive et al., 1998 
for 1E1724-3045), or alternatively by broken power laws 



(Wijnands and Van dcr Klis, 1999a, for example). Super- 
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Source 


Type" 


Flat topped noise and QPO 
References 13 


Red 


noise and QPO 
References 13 


4U 0614+09 


A 


"V 

i 


1 , 2 


Y 


44 c , 45 


4U 1608-52 


A 


Y 


<J A K 


_ 




4U 1702-42 


A 






Y 


14° 


4U 1705-44 


A 


"V 

i 


6, 5, 7 


Y 


5 C 


IE 1724-30 


A 


v 


8, 9 


_ 




4U 1728-34 


A 


v 

1 


10, 1 1 


_ 




KS 1731-260 


A 






Y 


15, 9 


4U 1735-44 


A 


Y 


9 


Y 


16 


SLX 1735-269 


A 


v 

I 


12, 9 


. 




4U 1744-26 


A 






Y 


17, 18 


XTE J 1806-246 


d 






Y 


19, 20 


SAX J1808. 4-3658 


P 


Y 


13 






GX 13+1 


A 






Y 


21 


4U 1812-12 


A 


Y 


2 






4U 1820-30 


A 






Y 


22 


GS 1826-238 


A 


Y 


f) 






4U 1915-05 


A 






Y 


46 


XTE J2123-058 


A 






Y 


23, 24 


Sco X-l 


Z 


Y 


25, 2 


Y 


25, 26, 27 


GX 340+0 


z 


Y 


2 


Y 


28 


GX 349+2 


z 






Y 


29 


GX 5-1 


z 


Y 


2 


Y 


30, 31, 32 


GX 17+2 


z 


Y 


2, 33 


Y 


33, 34, 35 


Cyg X-2 


z 


Y 


2, 36, 37, 38 


Y 


36, 37, 38, 39 


Cir X-l 


z e 


Y 


40, 41, 42, 43 


Y 


40, 42, 43 



a A, atoll source; Z, Z sou rce; P, millisecond X-ray pulsar. 
See Van Paradijs (1995) for complementary references up to 1995. 

c After the power spectra published in the literature. 

d Both atoll and Z classifications have been proposed ( Wijnands and Van der Klis, 1999b| ; Revnivtsev et al., 1999) 
c Shirey et al. ( [199$ }t999| ) 
References 



JL993) 



(4) Yu et al. (1997); 
(8) Olive et al. (1998); 



(1) Mndez et al. (|1997| ; (2) Wijnands and Van der klis (|1999a| ); (3) Yoshida et al. 
(5) Berger and V an de r Klis ( [1998| ); (6) Hasing er and Van der Klis ( |1989| ); (7) Ford et a l. ( |l998a| ) 

( 9) Ba rret et al. @000| ); (10) Strohmayer et al fll996| ); (11) Ford and V an de r Klis ( [l998| ); (12) Wijnands and V an de TKlis 
( |l999d) ; (13) Wij nands and Van der Klis ( [l998|); (14) Markwardt et a l. (|l99S| ); (15) Wijnands and Van der Klis (|l997 | ) ; (16) 
Wijnands et al. ( |l998c ); (17) Lewin et al. ( jl987|) ; (18) Strohmayer (|l99^~(19 ) Wijnands and Van der Klis (|l999b[ ); (20) 
Rev nivts ev et al. ( |1999 ); (21) Hom an et al. ( |l998j ); (22) Wijnands et al ( 1998| ); (23) Homan et al. ( |1999|); (2 4) Tomsick et 
( |l999|); ( 25) Van der Klis et al. ( |l997|) ; (26 TVan der Klis et al. ( |l99rj ); (27T Dieters and Van der K UsTPOOOI ) ; (28) Jonker 
(29) Kuulkers and Va n der Klis (|l998| ); (30) Kamadoet al. ( |l997| ); (31) Wijnands et al ( |l998b[ ); (32) Vaughan 
(33) Kuulkers et aL (^99^); (34) Wijnan ds et al. ( |l996| ); (35) Wijnand s et al. ( |l997a| ); (36) Wij nands et al. 
( |l997b| ); (37 ) Kuu lkers et al. ( |l999| ); (38) Focke ([l99(|; (39) Wijn ands et al. (|l998a|); (4 0) Oostcrb roek e t al. ( |l995| ); (41) 
Shirey et al. (|l996|); (42) Shirey et al. (|l99c|); (43) Shirey et al. (|l999|); (44) Ford et al. (|1997|); (45) Ford (|l997|); (46) This paper. 



Table 1. Sources for which a QPO (or a QPO-like feature or peaked HFN) has been reported in addition to either 
flat topped noise or red noise in their low frequency (~ 10~ 3 — 100 Hz) PDS. 



posed to the noise, a QPO-like feature is generally seen 
around - 0.5 - 70 Hz (see Table for a review). In this 
state, the root mean square (RMS) values of the fluc- 



tuations (to the mean flux) reach 50 % ( Van der Klis 
199|). On the other hand, on the banana branch, the 
shape of the PDS is roughly a power law with an index 
of -1.0; the s o-called very low frequency noise (VLFN) 
or red noise ( Hasinger and Van der Klis, 1989 ). In this 
state, the RMS is around 10 % or less. In the lower part 
of the banana branch, a broad structure, often refered 
to as "peaked HFN" , or a QPO is seen at frequencies in 
the range 1-70 Hz (see Table [j]). A correlation between 
the centroid frequency of the bump, and the frequency of 
the upper kHz QPO has been found in several LMXBs 
(e.g. KS1731-260), and interpreted in the framework of 



a Lcnsc-Thirring effect (Stella and Vietri, 1998; Stella 
et al., 1999j) or a two-oscillator model (Titarchuk and 
Oslerovich, 1999; Osherovich and Titarchuk, 1999a). 



1.2. 4U1915-05 

In the context of these results, we have analyzed our 
RXTE observations of the type I X-ray burster 4U1915- 
05. This source is known as the first dipping X-ray source 
discovered, providing the first reliable evidence for the bi- 

" White 



nary nature of such sources (Walter et al., 1982 
|and Swank, 1982| ; |5wank et al., 1984 ; [Bmale et al. 



1988) 



The period of the dips in X-rays is 50 minutes, the 



Walter et al., 1982 



ig sources ( 


White and Swank, 1982; 


Yoshida et al., 1995 




Chou et al.. 



1999 ). It is also o ptically identified with a 21st V mag- 
nitude blue star ( |Grindlay et al 



1 



. A modulation 
is seen in the optical but at a period ~ 1 % larger than 
the X-ray dip period ( Callanan et al., 1995 ). The discrep- 
ancy between the two periods is confirmed by Chou et al. 
(1999) who combined historical and recent data (includ- 
ing the RXTE data used in this paper) . This discrepancy 
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Name 


Start 




End 




P 




M 


Ttotal 


B 


D 


T 


R 


02/10 (1) 


Fob. 10 


00:14:04 


Feb. 10 


05:09:13 


5 




G 


6750 





3 


4112 


73.0 


03/13 (1) 


March 13 


22:41:24 


March 14 


00:54:13 


5 




G 


3316 








2880 


97.1 


05/05 


May 5 


05:59:39 


May 5 


09:46:13 


3, 


4 


G 


4039 








3904 


68.7 


05/06 


May 5 


22:16:00 


May 6 


03:37:13 


4 




G, CB 


6912 


1 





6592 


62.1 


05/14 


May 14 


09:22:05 


May 14 


12:55:13 


3 




E 


7599 





3 


4432 


29.5 


05/15 


May 15 


12:32:06 


May 15 


15:11:13 


3 




E 


6622 








6592 


41.5 


05/16 


May 16 


12:56:45 


May 16 


16:31:13 


3 




E 


7955 





4 


5936 


42.6 


05/17 


May 17 


08:20:08 


May 17 


11:49:13 


3 




G 


7640 





4 


2864 


35.4 


05/18 


May 18 


09:53:47 


May 18 


13:20:13 


3 




E 


7495 





4 


3168 


21.7 


05/19 


May 19 


10:06:38 


May 19 


13:37:13 


3 




E 


7389 





2 


2560 


22.6 


05/20 


May 20 


06:50:10 


May 20 


10:16:13 


3 




E 


7496 





2 


3904 


25.2 


05/21 


May 21 


07:47:52 


May 21 


11:22:13 


3 




E 


7949 





4 


2400 


23.3 


05/22 


May 22 


12:37:55 


May 22 


16:26:13 


3 




E 


8624 





6 


3584 


27.6 


05/23" 


May 23 


07:50:56 


May 23 


11:20:13 


4. 


5 


E 


7723 





8 


4128 


39.8 


06/01 


June 1 


17:36:53 


June 1 


21:36:13 


5 




E, G, CB 


9336 


1 


4 


4544 


46.3 


07/15 


July 15 


11:49:45 


July 15 


16:24:13 


5 




E 


9786 





6 


5168 


33.9 


08/16 


Aug. 16 


10:37:49 


Aug. 16 


14:55:13 


4. 


5 


E 


9983 


1 


5 


6384 


28.1 


09/06 


Sept. 6 


15:57:13 


Sept. 6 


21:26:13 


5 




E 


9344 





1 


6032 


39.5 


10/29 


Oct. 29 


06:53:21 


Oct. 29 


11:53:13 


5 




E 


9464 


1 


2 


5120 


35.9 



Table 2. Observation log. We list the name assigned to the observation, its start and stop times (day and UT hour), 
the number of PCUs working (P), the high time resolution data modes M available (G, E and CB designate the 
E_125us_64M_0_ls, GoodXenon_16s and CB_2ms_64M_0_249 (Binned Burst Catcher) configurations respectively), the 
total exposure time (T to tai) in seconds after screening (on elevation and offset), the number of bursts (B) and dips (D) 
present, the exposure time in seconds of the persistent light curve after filtering out the bursting and dipping parts 
(T) and the 2-20 keV background subtracted count rate in units of cts s _1 PCU -1 (R) of the persistent emission. 
W Gain epoch 1 data (the other observations are gain epoch 3). ° Two small observations were carried out on May 
23rd, we merged them together for clarity. 



has led to several interpretations among which the SU 
UMaQ model is prefered, although the hierarchical triple 
system model proposed by Grindlay et al. ( [1988 ) is still 
possible ( Chou et al., 1999| ). 



Gamma- Ray Observatory had failed to detect the burster 



To date, the only timing study of the non-dipping 
and non-bursting emission of 4U1915-05 was carried out 
with the Ginga data, which represented the most exten- 
sive X-ray observations of the source before the launch of 
BeppoSAX and RXTE. The 1-37 keV luminosity varied 
between 5 x 10 36 and 1.2 x 10 37 ergs s" 1 , in the range of lu- 
minosities previously observed for the source. A V-shape 
was observed in th e hardness (7-1 8/1-7 keV) intensity (1- 
18 keV) diagram ( |Yoshida, 199^) . The 3 x 1CT 3 - 1 Hz 
PDS were fitted with a power law of index -0.6 for the 
left branch of the V data an d of i ndex roughly -1.1 for the 
right branch data. Yoshida ( 1992) sug gested that 4U1915- 
05 was an atoll source. Yoshida ( 1992] ) interpreted the left 
branch of the hardness intensity diagram as the island 
state and the right branch as the banana state, although 
the two branches overlapped in the color-color diagram. 

From an X-ray burst showing photospheric radius 
expansion in the GINGA data, the distance to the 



source was estimated as 9.3 kpc (Yoshida, 1992). 4U1915- 
05 was observed at energies up to 100 keV by BeppoSAX 
(Church et al., 1998). This was the first observation of a 
hard X-ray tail in the spectrum of 4U1915-05. Indeed, 
previous pointed and monitoring observations at high 
energies with OSSE and BATSE aboard the Compton 



1 SU UMa stars are cataclysmic variables, which during 
their long ('super') outbursts show humps which have a period 
which is a few % longer than their orbital period. 



( [Barret et al., 1996| ; |Bloser et al., 19960 . 4U1915-05 is the 
first dipping source to belong to the growing group of 
LMXBs detected up to - 100 keV ([Barret et al., 2000|). 



In this paper we report on the timing analysis of the 
X-ray emission which led to the discovery of HFQPOs 
and LFQPOs in 4U1915-05. In the next section, we de- 
scribe the observations and the data analysis. Then, we 
present the results. Finally, we compare 4U1915-05 to 
similar sources and discuss the results in the framework 
of currently proposed models of HF and LFQPOs. 

2. Observations and data analysis 

2.1. The RXTE/PCA observations 

4U1915-05 was observed by RXTE on 19 occasions in 
1996: on February 10th, March 13th, May 5-6th, May 
14th to 23rd, June 1st, July 15th, August 16th, Septem- 
ber 6th and October 29th. The log of the observations 
is presented in Table |[ An observation has an exposure 
time of roughly 8000 seconds and is provided as typi- 
cally three data files with gaps between files. These files 
are refered as segments of observation in our analysis. 
We report on the data collected in the 2-60 keV range 
with the PCA. It consists of 5 nearly identical large area 
Proportional Counter Units (PCU to 4) corresponding 
to a total collecting area of ~ 6500 cm 2 ( Jahoda et al 



1996 ). For safety reasons, PCUs are switched on and off 
in the course of an observation. During our observations, 
the data were mainly obtained in three different PCU 
configurations: all PCUs on, first four on and first three 
on. Furthermore, our observations cover two PCA gain 
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100 150 200 

Time (Days after 1996/02/24) 



Fig. 1. Light curve of 4U1915-05 as observed by the 
RXTE/ ASM. The binning time is one day. The times 
of the pointed observations are indicated by arrows. 



epochs: the February 10th and March 13th observations 
belong to the first one whereas the following observations 
belong to the third one. Different response matrices were 
used for each set. More than 138 ks of good data are avail- 
able after recommended screening (elevation angle above 
the Earth's limb greater than 10 degrees and pointing 
offset angle less than 0.02 degree). 

The all-sky monitor (ASM) aboard RXTE consists of 
three coded-aperture cameras mounted on a motorized 
rotation drive, allowing to view different regio ns of the 
sky during a satellite orbit ( Levine et al., 1996 ) . The in- 
strument is sensitive in the energy range ~ 2-12 keV. Fig. 
| is the ASM light curve of 4U1915-05 during 1996. The 
binning time is one day. The times of the PCA obser- 
vations are flagged with arrows. Four X-ray bursts were 
recorded during the PCA observations: on May 5th, June 
1st, August 16th and October 29th. With the exception 
of the first one, they began during a dip. A total of 58 pri- 
mary or secondary dips (or parts of dips) were observed, 
lasting typically 400 seconds. Dips are thought to be due 
to occultation of the internal emitting region by vertical 



structure at the outer edge of the accretion disk (Walter 
E 1982fc |White and Swank, 1982]) 



et 



Primary dips oc- 
cur at intervals consistent with an orbital period of ~ 50 
minutes. Secondary dips are "anomalous" dips, narrower 
than the primary ones. They occur irregularly at ~ 180° 
out of phase with the primary dips ( Chou et al., 1999| ). 



Dipping is almost 100 % in the 3-5 keV energy range, but 
only ~ 65 % in the 5-30 keV range. 

2.2. Data modes 

The Standard 2 mode is available for each observation. 
This mode provides counts integrated during 16 seconds 



in 128 energy channels covering the 2-60 keV energy 
range. The high time resolution data were provided in two 
different modes, depending on the observations (see Ta- 
ble |): the E_125s_64M_0_ls mode with 122 fis resolution 
or the GoodXenon_16s mode with ~ 1 /j,s resolution. The 
bursts that occured on October 29th and August 16th 
were recorded in the 122 [is resolution mode. The bursts 
that occured on May 5th and June 1st were recorded 
during 3.75 seconds in a burst catcher mode with 2 ms 
resolution. 



2.3. Data analysis 

2.3.1. Light curves, color-color and hardness-intensity 
diagrams 

For the light curves, color-color and hardness-intensity 
diagrams, we have used the Standard 2 mode data. We 
made light curves in adjacent energy bands. The bursts 
and dips were then filtered out after visual screening of 
the light curves to generate the color-color and hardness- 
intensity diagrams for the persistent emission. The soft 
color was calculated as the 3-5 keV/ 1.7-3 keV count rate 
ratio, and the hard color as 10-30 keV/5-10 keV. 

As all the PCUs have different energy responses, the 
absence of a PCU can mimic spectral variations. Thus, 
we used the data collected with the first three PCUs since 
they were continuously operating during all our observa- 
tions. 

Differences between response matrix of gain epoch 1 
and 3 introduce differences in the boundaries of the en- 
ergy channels used. This instrumental effect introduces 
shifts between the colors or count rates obtained from dif- 
ferent epochs. We quantified these shifts using the Crab 
nebula which is a steady source. We find that the cor- 
rection factors to apply to gain epoch 1 data in order to 
make them match the epoch 3 region are +10.1 %, -1.5 %, 
-1.1 % and +2.8 % for the count rates in the energy bands 
1.7-3, 3-5, 5-10 and 10-30 keV respectively. The correc- 
tion factors are -10.5 % and +4.0 % for the soft and hard 
colors respectively. We applied these correction factors to 
4U1915-05 data. 



2.3.2. Power density spectra 

To investigate the variability of 4U1915-05, we computed 
PDS using the FTOOLS powspec. We treated separately 
the data obtained from the different high time resolution 
modes. To reach timescales ranging from roughly 10 ms to 
4 minutes, the data were rebinned to a ~ 4 ms resolution. 
Then, each continuous set was divided into segments of 
65536 bins. A fast Fourier transform (FFT) of each seg- 
ment was computed yielding a low frequency PDS in the 
range - 4 x 10~ 3 - 128 Hz. 

To reach timescales ranging from roughly 0.25 ms to 1 
second, the data were rebinned to a ~ 122 fj,s resolution. 
Each continuous set was divided into segments of 8192 
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bins. A FFT of each segment was computed leading to 
high frequency PDS in the range 1-4096 Hz. Then, the 
low or high frequency PDS obtained from a given set of 
observations, a given observation or a given segment of 
observation were averaged together to obtain a final PDS 
representative of the set, observation or segment. 

The PDS were normalized according to Leahy et al. 



quency bins, the error on the binned points becomes 



E F 



VMWN 



and the SNR can be computed as: 



(1983) so that the average power expected from a Poisson 
distribution is 2. We have checked in the high frequency 
PDS that the average white noise level above 1500 Hz 
was indeed 2, indicating that deadtime corrections were 
not needed. 

The error on each PDS bin has been set to Ep = 
f™ where M is the number of raw PDS averaged to- 
gether, W the number of raw frequency bins averaged to- 
gether (Van dcr Klis, 1989). P m is set to the Fourier power 
at the point considered for low frequency PDS (where the 
power can differ significantly from the white noise level 
2 because of the source aperiodic variability) and to 2 
for the high frequency PDS. Both logarithmic and linear 
rcbinnings were applied. 

Different components were summed to fit the result- 
ing PDS: a constant for the white noise level (C), a power 
law for the VLFN (PL), a cutoff power law for the HFN 
(CPL). The functional shapes used for the latter com- 
ponents are kv a and Ai/ a exp(— f/f cutoff) respectively, 
where A is the normalization, a the power law index, 
v the frequency and Cutoff the cutoff frequency. C is a 
free parameter for the high frequency PDS and is set to 2 
for the low frequency PDS. We used Gaussians (G) to de- 
scribe the QPO features. Note that Lorentzians are also 



currently used (Van der Klis, 1995). 

All the fits were performed using a \ 2 minimization 
technique. The errors on the analytic model parameters 
correspond to a \ 2 variation of 2.7, equivalent to the 90 % 
confidence region for a single interesting parameter. 

To compute the RMS, the PD S wer e normalized ac- 
cording to Bclloni and Hasinger (199C) taking into ac- 
count the total and background count rates. The RMS 
was then obtained by integrating the normalized PDS in 
the appropriate frequency range. 

We have also computed PDS of data from the bursts. 
For the bursts recorded in the 122 fis resolution mode, we 
have calculated an FFT power spectrum using 4096 bins 
lasting 244 fj,s, so that one PDS corresponds to a segment 
of the burst of 1 sec duration and the Nyquist frequency 
is 2048 Hz. For the other bursts, the time resolution is 
only ~ 2 ms, yielding a Nyquist frequency of 256 Hz. 

2.3.3. Search technique for quasi-periodic oscillations 

To look for weak QPOs in the PDS (as in the case of 
4U1915-05), we developed an algorithm which computes 
for each PDS the width of the window which maximises 
the signal to noise ratio (SNR) within the window. As- 
suming that a given window of width w contains N fre- 



SNR{in cr) 



N S Pj P 



ref 



E P 



(1) 



where Pj are the N PDS points contained in the win- 
dow of width w and P r ef is the reference power. For 
the low frequency PDS, given the weakness of the VLFN 
component above ~ 1 Hz, we looked for excesses above 
P re f=2, the theoretical level expected from Poissonian 
noise. PDS have been processed through the algorithm 
and searched for excesses of widths from 1 to 50 Hz in 
the frequency range 1-100 Hz. For the high frequency 
PDS, we looked for excesses of widths 5-200 Hz in the 
frequency range 100-1500 Hz. We used the noise level av- 
eraged above 1500 Hz (~ 2.0) as the reference P re f for 
the search of excesses in the PDS. In both cases, we first 
selected the positive detections above 3.5 cr. A new signif- 
icance level inferred from the fit was then attributed to 
each excess in order to take into account the excess shape. 
We report here on the LFQPOs and HFQPOs which have 
a significance level greater than 3 cr. This significance level 
is reported in Table |] with the properties of the QPOs 
detected in 4U1915-05. We do mention either the obser- 
vation or the segments, so that no overlapping is possible. 

This algorithm does not take into account the num- 
ber of trials in determining the significance level. To de- 
termine the r eal si gnificance of the detections, we follow 
Van der Klis ( 1989 ) and define the (1 — e) confidence level 
P detect as the power level that has only the small probabil- 
ity e/N tr j a ; to be exceeded by a noise power. N tria ; is the 
number of different power values that one wishes to com- 
pare with P detect- It may be the total number of bins in 
the PDS, or less if only a given frequency range is looked 
for excesses. The number of different PDS looked at may 
also be included in the number of trials but we did not 
take it into consideration. In our case, the high frequency 
PDS are scanned over a frequency range of 1400 Hz. The 
exposure time of a raw PDS is 1 second (corresponding 
to a frequency resolution of 1 Hz). The typical exposure 
time for a segment of observation is 3000 seconds so that 
M=3000 raw PDS are averaged together to create the 
PDS representative of the segment. Assuming that the 
algorithm finds a maximum excess of width 20 Hz (cor- 
responding to a final rebinning factor NW of 20) at 5 cr, 
then N tr j O j=70 and the actual significance level of the 
signal is 99.998 %. In the worst case of an excess of width 
5 Hz (ie the minimum width considered) found at 3.5 cr 
(ie the minimum level obtained from the algorithm con- 
sidered), N tr i a ;=280 and the actual significance level is 
92.7 %. The low frequency PDS are scanned over a fre- 
quency range of 100 Hz. The exposure time of a raw PDS 
is ~ 256 seconds (corresponding to a frequency resolution 
of 4 x 10~ 3 Hz). A PDS representative of a 3000 seconds 
segment is obtained from M=12 raw PDS. In the worst 
case of an excess of width 1 Hz (e.g NW=256) at 3.5 cr, 
then N tria / = 100 and the true significance level is 97.0 %. 
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3. Results 

3.1. Light curves, color- color and hardness-intensity 
diagrams 

Fig. H shows the global evolution of the persistent emis- 
sion during our observations. Each point gives the 2- 
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Fig. 2. Evolution of the persistent emission during the 
observations. The 2-20 kcV background subtracted count 
rates are averaged over each observation. The correspond- 
ing exposure time is given in Table |^. The symbols labeled 
with the date of the observation will be used throughout 
the remainder of this paper. 



20 keV background subtracted count rate integrated over 
a single observation. Each plotting symbol is associated 
with an observation (see the right corner of Fig. |2|). As 
could have been anticipated from Fig. |l|, our observations 
sampled various intensity states of the source: the PCA 
count rate changes by a factor of ~ 5 from roughly 20 
to 100 cts s _1 PCU -1 . We also generated light curves 
in different energy bands between 2 and 30 keV for each 
observation. No substantial variation is noticeable in the 
persistent emission over timescales of a few hours. 

Fig. [3] shows the color-color diagram. The points seem 
to form a "Banana" shape than can be approximated by 
a spline. 
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Fig. 3. Color-Color diagram. The binning time is 512 sec- 
onds and a typical error bar is shown. The thick line is 
the spline used to approximate the position on the color- 
color diagram. This position is quantified by the curvilin- 
ear coordinate S a chosen arbitrary. The positions S a = l 
and S a =2 are indicated by arrows. 
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Fig. 4. Hardness-intensity diagrams for the soft color 
(top panel) and the hard color (bottom panel). In the 
main panels, the binning time is 512 seconds and a typ- 
ical error bar is shown. In the inset panels, each point is 
the average over an observation. The asterisks and the 
large crosses indicate observations of the LS and HS re- 
spectively 
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Frequency (Hz) Frequencj' (Hz) 

Fig. 5. Average PDS in the LS (left panel) and HS (right panel). The first one is obtained by averaging the PDS of 
the LS observations available in the 122 fj,s resolution mode: all the LS observations except a segment of the June 1st 
observation and the May 17th observation. The second PDS is obtained by averaging the PDS of the HS observations: 
February 10th, March 13th, May 5th and 6th (all obtained in 1 fis resolution mode). The total exposure time of the 
first PDS is more than 15 hours whereas it is only ~ 3 hours for the second one. The inset panel is a zoom of the 
PDS above 1 Hz. The solid lines are the best fits to the data. 
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Table 3. Best fit parameters obtained for the state averag ed PDS shown in Fig. |. The RMS is given in the frequency 
range 4 x KT 3 - 1 Hz for the VLFN, 1-128 Hz for the HFN and 0.1-0.3 Hz for the QPO. 



Fig. [| shows the soft (top panel) and hard (bottom 
panel) colors as a function of the count rate. In the main 
panels, the binning time is 512 seconds while in the inset 
panels, each point is the average over an observation. 

Looking at the behaviour of the source within an ob- 
servation, we see, from the vertical tracks drawn by the 
points with a same symbol (main panels) , that the colors 
can vary a lot while the intensity remains nearly constant, 
especially the hard color at low count rates. 

Looking at the global behaviour of the source along 
the observations (inset panels), two groups correspond- 
ing to two different spectral states may be distingished: 
we refer to the observations with the 5-30 keV intensity 
below and above ~ 30 cts s _1 PCU -1 as being in the low 
and high states (LS, HS) respectively. In the LS, we see 
a trend for the mean hard color to decrease as the inten- 
sity increases (Fig. [| bottom), while the soft color remains 
nearly constant (Fig. [|top), so that, on average, the emis- 
sion softens with the intensity. On the contrary, in the HS, 
the hard color is nearly constant (Fig. |] bottom), while 
the soft color increases with the intensity (Fig. ^ top), 
so that, on average, the emission hardens. 4U1915-05 is 



in a low intensity/hard spectrum regime in the LS and 
in a high intensity/soft spectrum regime in the HS. The 
spectral modeling of the pe rsistent emission confirms this 
picture ( B loser ct al., 200C ) . Assuming a source distance 
of 9.3 kpc, we derived a 2-20 keV luminosity of 1.4 x 10 37 
ergs s _1 on March 13th (when 4U1915-05 is the brighest 
and in the HS) and of 3.2 x 10 36 ergs s" 1 on May 18th 
(when 4U1915-05 is the faintest and in the LS). 



3.2. Variability in the 4 X 10 -3 — 128 Hz frequency range 

3.2.1. Average power density spectra 

In order to better characterize the source, we looked 
for possible correlations between the aperiodic variability 
and the spectral behaviour. We thus averaged together 
low frequency PDS of the LS on one hand (Fig. [3] left), 
and PDS of the HS on the other hand (Fig. || right), 
provided that they were obtained from the same timing 
mode data. 

The average PDS of Fig. || (left) has an integrated 
RMS over the entire ~ 4 x 10~ 3 — 128 Hz frequency range 
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of 18.8±0.6 %. The PDS shows VLFN up to roughly 1 Hz 
and HFN around 20 Hz. The model C+PL+CPL provides 
a better fit than a simple PL (x 2 /d.o.f ~ 166/146 and 
324/f49 respectively). Looking at Fig. || (right), we note 
an excess power around 0.2 Hz. We have fitted that ex- 
cess with a Gaussian and assessed its significance through 
the standard F-test. We define P as the probability of 
rejecting the hypothesis that the fit is better with the 
additional Gaussian component. In the present case, we 
have P = 0.09, which means that the QPO is significant 
at a level of ~ 9f %. With the 0.2 Hz Gaussian added, 
the results of the best fit are given in Table ||. 
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Frequency (Hz) 
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Fig. 6. PDS showing LFQPOs in the persistent emission 
of two segments of observations on October 29th (top 
panel) and June 1st (bottom panel). The mean count 
rate is 35 cts s" 1 PCU" 1 for the first PDS and 41 cts s" 1 
PCU -1 for the second (no background subtraction). The 
exposure time of both PDS is 1800 seconds. The inset 
panel is a zoom of the PDS above 1 Hz. The solid lines 
are the best fits to the data. 



We have searched for the same signal in different en- 
ergy ranges. As an example, in the 2-30 keV band, the 
significance of the 0.2 Hz signal reaches 93 %. The best 
fit parameters are given in Table ^. In addition to the 
state averaged PDS, we computed PDS for each observ- 
ing day. In some of these PDS, power excesses are visible 
around 0.2 Hz. However, they are statisticaly marginal 
with significance levels less than 90 %. 

In the HS, the source displays very low variability. 
The total RMS in the 5-30 keV energy range is 5.3 ± 
0.2 %. At low frequencies, the PDS is still dominated by 
VLFN up to roughly 0.1 Hz. But above this frequency, 
the PDS is consistent with having no power. No QPO 
nor HFN component are present in the PDS. The best fit 
parameters are given in Table ||. 

Looking at each individual PDS, we found that what 
we called previously the HFN component in the state 
averaged PDS could be in fact resolved into broad QPO- 
like features of varying centroid frequencies. In the fol- 
lowing, we call these features LFQPOs. Fig. ^ shows two 
examples of daily-averaged PDS where such a LFQPO is 
clearly detected. Detailed properties of the LFQPOs are 
given in the next section. To conclude, 4U1915-05 dis- 
plays aperiodic variability below ~ 1 Hz (the VLFN) and 
quasi-periodic variability or no variability above 1 Hz. 

Fig. displays the RMS as a function of the count 
rate for each observation. The RMS is given in the 3 x 
10~ 3 — 1 Hz frequency range where the VLFN dominates, 
and between 1-100 Hz where the LFQPO is detected. The 
RMS of the VLFN component ranges between ~ 3 and 
6 %, whereas the RMS above 1 Hz ranges between ~ 5 
and 25 %. In the LS, the RMS is high and decreases as the 
average count rate per observation increases. On the other 
hand, in the HS, the RMS, after reaching a low threshold, 
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Fig. 7. RMS in the 4 x 10~ 3 — 1 Hz range (small symbols) 
and in the 1-128 Hz range (big symbols) versus count rate 
(no background subtraction) . 
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does not significantly vary as the count rate increases. 
Thus, the timing behaviour of 4U1915-05 is different in 
its two states. 



3.2.2. Low frequency quasi-periodic oscillations 

To further investigate the properties of the LFQPOs de- 
tected between a few Hz and 100 Hz in 4U1915-05 (see 
Fig. and given their relative weakness in the persis- 
tent emission, we searched for them in the persistent and 
dipping emission combined rather than in the persistent 
emission alone. Although we note that we may be look- 
ing at different signals when combining the persistent 
and dipping emissions, we aimed to get the best signal 
to noise ratio. The study of the LFQPOs is carried out 
in the 5-30 keV energy range in order to make easier 
the comparison with the HFQPOs detected in that en- 
ergy range with the best signal to noise ratio (see next 
section). The dipping is large but not complete in that 
energy range. So, adding the dipping parts, which rep- 
resent roughly 40 % of the data set, actually increases 
the signal to noise ratio. This is not negligible because 
4U1915-05 is relatively faint in the PCA (2-20 keV back- 
ground subracted count rate around 35 cts s _1 PCU -1 , 
to be compared with ~ 30000 cts s" 1 PCU" 1 for Sco X-l) 
and observed with only 3 PCUs for more than the half of 
our observation (see Table ||). Furthermore, the ingress 
or egress of a dip lasts more than 1 second (typically 40 
s). So, above 1 Hz where the LFQPOs are detected, the 
PDS of the persistent and dipping emission is not contam- 
inated by the dip phenomenon itself. PDS were computed 
in the range ~4x 10~ 3 — 128 Hz the same way as for the 
persistent emission alone. We averaged together PDS of 
the same observation or the same segment. 

The LFQPOs frequency, FWHM, RMS and signifi- 
cance level are listed in Table ||. Their frequency ranges 
between 5 and 80 Hz. Their FWHM is around 16 Hz 
but can vary from 0.5 to 45 Hz, so that their coherence 
(Q=y/FWHM) changes roughly from 0.5 to 26 with a 
mean value of 4. The designation of QPO is then some- 
times improper since this is usually used for signals with 
Q values greater than 2 (Van der Klis, 1995). One of the 



LF QPOs is shown in Fig. |3 . In the second segment of 
September 6th, we found an excess in the PDS between 
1 and ~ 50 Hz (~ 6.3 a). A visual inspection of the PDS 
indicates that this excess is flat, yet above the noise level. 
For this reason, we were unable to fit it with a Gaussian, 
and therefore we do not consider this excess as a LFQPO, 
as defined above. Interestingly enough, however we note 
that the centroid of the excess (~ 30 Hz) would fall on 
the correlation between LFQPO frequency and the upper 
HFQPO frequency, as shown in Fig. [jj] (see section 4.4). 

We looked for correlations between LFQPOs and 
spectral states. The bottom-left panel of Fig. || shows 
the LFQPO frequency versus the 5-30 keV count rate. 
All but one of the LFQPOs are detected in observations 
of the LS. For 41X1915-05, the LFQPO frequency seems 
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Fig. 8. PDS showing the LFQPO in a segment of the Oc- 
tober 29th observation (persistent and dipping emissions 
combined). The solid line is the best fit of the QPO with 
a Gaussian. 



to roughly correlate positively with the intensity, but the 
scatter is substantial. 

The bottom-right panel of Fig. ^| shows the LFQPO 
RMS as a function of count rate. The RMS ranges from 5 
to 20 %. For the HS, we have derived upper limits ranging 
between 3 and 10 % (see Table |j). No simple relation 
seems to connect the two quantities. 

Recent works have shown that some timing features 
of similar systems were better correlated with the posi- 
tion on the color-color diagram than with the count rate 
flMendez et al., 199S| ; pVlendez and Van der Klis, 1999Q . 
We thus examined the QPO properties as a function of 
the position on the color-color diagram. 

Color-color diagrams are sensitive to changes of a few 
percent in flux ratios between energy bands. Observations 
allowed to show that during these subtle changes, the 
X-ray intensity was not a measure of the accretion rate 
since both positive and negative correlations between the 



two quantities could occ ur (Hasinger and Van der Klis 



1989J [Van dor Klis, 1994[) . The position on the color-color 



diagram is a better indicator of the spectral state and 
hence the accretion rate than the count rate. 

The direction of evolution of the accretion rate on the 
color-color diagrams was inferred from X-ray luminosities 
and c onfirmed by multiwavclcngth and bursts obse rva- 
tions ( Hasinger et al., 1990 ; Van der Klis et al., 1990 ). In 
the atoll sources, the accretion rate is believed to increase 
from the island to the lower and then to upper banana 
states. For our analysis, we refer to the color-color dia- 
gram shown in Fig. pi The spline shows the approximated 
track followed by 4U1915-05 during the observations. The 
direction of evolution of the accretion rate along the track 
is derived from the X-ray luminosity which globaly in- 
creases from the extremity located at the top of the dia- 
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Table 4. Properties of the QPOs detected above 3 a (see text) in 4U1915-05 in the 5-30 keV energy range. We 
list the observation date (month/day) followed by the segment number or by an asterisk if all segments have been 
summed, the frequency, FWHM, RMS and significance of the LFQPOs (1-100 Hz) and HFQPOs (100-1400 Hz). 
The last column is the separation between the twin HFQPOs. All the values are inferred from Gaussians used to fit 
the QPOs. The reduced x 2 are close to 1. The error bars are given for the 90 % confidence level (Ax 2 =2.7 for the 
variation of a single parameter). When no signal is present above our confidence threshold, we give a 3 a upper limit 
on the RMS of a signal of FWHM 10 Hz for the LFQPOs and 50 Hz for the HFQPOs. The significances quoted here 
are for a single trial only. A more detailed assessment of the significance of the detections taking into account the 



number of trials may be found in the text (see section 2.3.3). a flags detections at a significance level between 3 and 
4 a. The other signals have significance level greater than 4 a. b The excess could not be fitted as a Gaussian (see 
text). Approximated values of the frequency, FWHM and RMS are given. 



gram to the extremity located at the bottom (see Bloser 
ct al. (2000) for details about spectral fitting). 



Following Mndez et al. ( |1999| ), we called S a the 
parametrization along this spline. S a is a measure of the 



position of the source on the diagram. S a was choosen to 
increase from top to bottom so that it is representative 
of the inferred evolution of the accretion rate. We deter- 
mined the value of S a for each segment of observation by 
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Fig. 9. Left: LFQPO (bottom panel) and HFQPO (top panel) frequency as a function of the 5-30 keV count rate (no 
background subtraction). Right: LFQPO (bottom panel) and HFQPO (top panel) RMS as a function of the 5-30 keV 
count rate. Small symbols are used for detections between 3 and 4 a, big symbols for detections above 4 a. 



projecting the corresponding point of the color-color dia- 
gram (inset panel of Fig. where one point corresponds 
to one segment) on the spline. This value of S Q could then 
also be attributed to the LFQPO detected in the given 
segment. For LFQPOs detected in PDS averaged from 
several segments of the same observation day, we used 
the mean of the S a values found for each segment. The 
bottom-left panel of Fig. [l(] shows the LFQPO frequency 
as a function of S a . A positive correlation is now clearly 
visible. In the bottom-right panel of Fig. we show that 
the LFQPO RMS anticorrelates with S a . In both plots, 
the scatter is now much smaller than in the previous plots 
where the count rate was involved. Note that we have also 
studied the evolution of the LFQPOs as a function of the 
hard color (10-30 keV/5-10 keV). We concluded that the 
hard color was a worse indicator of the timing behaviour 



than S a : the evolution of the soft color had to be taken 
into account. 



3. 3. Discovery of high frequency quasi-periodic 
oscillations 

Now, we report on the discovery of HFQPOs from 
4U1915-05. The analysis has been carried out in the 5 
to 30 keV energy range because in other sources, it is 
known that HFQPOs are more easily detected in this en- 
band, as in, e.g. KS1731-260 (Wijnands and Van 



ergy 



der Klis, 1997). The persistent and dipping emission were 
combined to increase the signal to noise ratio (see section 
l2j). 



The results are reported in Table [| (either the obser- 
vations or the segments). The HFQPOs frequency ranges 
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Fig. 10. Left: LFQPO (bottom panel) and HFQPO (top panel) frequency as a function of the position S a on the 
color-color diagram. The asterisks represent the frequency separation of the twin HFQPOs. The solid line is the 
average separation. Right: LFQPO (bottom panel) and HFQPO (top panel) RMS versus S a - Small symbols are used 
for detections between 3 and 4 a, big symbols for detections above 4 a. 



between - 200 and 1300 Hz. Their FWHM is around 
70 Hz and coherence around 50. Their RMS amplitude is 
high (~ 17 %). One of the strongest HFQPO signals is 
shown in Fig. [ll]. 

We detect twin HFQPOs in 5 segments or observa- 
tions (Table ||). For four of them, the frequency separation 
is consistent with being constant. The mean separation is 
348 ± 12 Hz. For the fifth, the separation is 290 ± 5 Hz, 
thus inconsistent with the previous value. We note that 
this pair of QPOs occur at 514 and 224 Hz whereas the 
four other twin peaks occur above 550 Hz. The twin HFQ- 
POs detected in a segment of the June 1st observation are 
shown in Fig. [l^. 

HFQPOs were not detected in observations with the 
largest count rates (February 2nd, March 13rd) with up- 



per limits of - 6 % on the RMS (Table |). HFQPOs were 
not detected either in the lowest count rate regime (May 
18th to May 22nd); this might be due to a lack of sen- 
sitivity. Indeed, only 3 PCUs were working during these 
observations (see Table |^) and the upper limits on the 
RMS are not constraining (~ 20 %, see Table ||). 

The top-left panel of Fig. | shows the HFQPO fre- 
quency as a function of count rate. In this diagram, we 
observe a branch between 700 and 1100 Hz where the fre- 
quency strictly correlates with the count rate. A second 
parallel branch seems to be drawn at lower frequencies. 
However, some of the HFQPOs of that branch do not 
have a simultaneous twin peak in the upper branch. In 
the lowest count rate regime, the frequency does not fol- 
low anymore a simple relation with the count rate. 
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Fig. 11. PDS for the September 6th observation. 
HFQPO is fitted with a Gaussian. 
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Fig. 12. The twin HFQPOs fitted with Gaussians de- 
tected in a segment of the June 1st observation. For dis- 
play purposes, two different frequency binnings have been 
used. 



The top-left panel of Fig. |10| shows the HFQPO fre- 
quency as a function of S a . The frequency versus S a re- 
lation appears much simpler than the frequency versus 
count rate relation. We can now easily identify the upper 
and lower HFQPOs. The frequency of the upper HFQPO 
is now well correlated with S a within its full range. The 
four twin peaks above 550 Hz are detected within a nar- 
row range of S Q around 3. Their frequency separation is 
shown with an asterisk on Fig. [Io|. The pair of QPOs at 
224 Hz and 514 Hz appears at the lowest S a and hence 
lowest inferred accretion rate. Concerning the HFQPO 
RMS, we show that, surprisingly, it is better correlated 



Fig. 13. LFQPO frequency as a function of the upper 
HFQPO frequency. We report signals detected simulta- 
neously above 4 a. The solid line is the powerlaw of 
index 1.6 fitting the data. The dashed line is the pre- 
cession frequency (Lense-Thirring and classical contribu- 
tions (Stella and Vietri 1998)) assuming standard param- 
eters for the neutron star (see text) and a spin frequency 
of 348 Hz. 



(negatively) with the count rate than with S Q (top-right 
panel of Fig. ^ and [l0| respectively). 

3.4- Correlation between low and high frequency 
quasi-periodic oscillations 

Fig. [l(] (left panel) shows the parallelism between HF and 
LF QPOs. Their frequency increases similarly with S Q . In 
Fig. |l3|, we show the relation between LFQPO and up- 
per HFQPO frequencies when both QPOs are detected 
simultaneously and above 4 a. The correlation between 
the two frequencies is obvious. It can be approximated 
with a power law of the form vlfqpo = av HFQPO- ^ ne 
first attempt of fitting the data provided a high reduced 
X 2 because of the point representing the July 15th obser- 
vation lying off from the fit. We thus made the fit taking 
out this point. We obtained a = (6 ± 10) x 10~ 4 and 
(3 = 1.6 ± 0.3 (reduced x 2 °f 1-1) as parameters of the 
power law drawn as the solid line on Fig. 131 



3. 5. Search for oscillations in bursts 

For the bursts recorded in the 122 fis resolution mode, we 
have computed thirty PDS of 1 second duration so that 
they covered the rising and decaying parts of the bursts. 
Each was individually searched for excess power. No such 
excesses were found either in individual PDS nor in their 
sum. We therefore conclude that no QPOs were found 
in those bursts. We have set upper limits on the RMS, 
assuming a 3 a signal of FWHM 1 Hz in the PDS. We 



14 



Boirin et al.: RXTE observation of 4U1915-05 



derive a 3 a upper limit of 3.2 % for the October 29th 
burst, and of 2.8 % for the August 16th burst. 

For the burst catcher mode, the situation is slightly 
different. The data cover only the first 3.75 seconds of 
the bursts, and the Nyquist frequency is 256 Hz which is 



below the typical frequency of a burst oscillation (Van der 
Klij, 1999| ) . Note however that coherent oscillations were 



found at 7.6 Hz in a burst from Aql X-l ( {Bchoclkopf and 



Kel [ey, 1991 ). For the burst that occured on May 5th 



no significant signals were found in PDS computed for 
different segments durations (ranging from 1.25 to 3.75 
seconds). We have derived a 3 a upper limit of 2.8 % 
(FWHM=1 Hz) for the 3.75 seconds burst duration. 

For the June 1st burst, we have also computed an 
FFT from the 3.75 seconds data set. When the PDS so 
computed is rebinned by a factor of 4 (Av — 1.1 Hz), two 
signals are detected by our algorithm; 117.5 Hz (5.4 a, 
RMS=3.5 ± 0.4 %), 225.7 Hz (5.0 a, RMS=2.8 ± 0.2 %). 
However, in the case of these PDS, the corresponding sig- 
nificance detection l evel (t aking into account the number 
of trials; see section 2.3.3 ) is only 85 %, which is insuffi- 
cient to claim a detection. 

4. Discussion 



is the typical aperiodic variability of atoll sources in their 
island state; 4U1915-05 displays VLFN in its both states. 
This behaviour was observed in the GINGA data as well. 
Unlike Yoshida ( 1992| ) who associated the two regimes 
with island and banana states, we believe that, during 
our observations, the source was not in its island state; it 
probably remained on the banana branch and hence never 
reached the state which may be at even lower luminos- 
ity. Such a state may have been observed by BeppoSAX 
( Church et al., 1998] ) since extrapolati ng the spectral pa- 
rameters reported by Church et al. (fL998) into the 2- 



50 keV energy range gives a luminosity of 4.5 x 10 ergs 
s _1 while the luminosity derived from the RXTE obser- 
vations in the same ener gy band ranges bet ween 5 x 10 36 



and 1.4 x 10 37 ergs s" 1 ([Bloser et al., 2000|) . 



The properties and presence of HFQPOs strongly de- 
pend on source states. For atoll sources, HFQPOs seem 
to occur at intermediate inferred accretion rates: they are 
generally not observed in extreme island or upper banana 
states corresponding respectively to the lowest and high- 
est accretion rates in a given source, as in e.g. 4U1608-52 



(Mendez et al., 199E). Twin simultaneous HFQPOs have 
now been detected in all atoll sources showing kHz vari- 



ability but Aql X-l ( |Van der Klis, 1999j ). They are always 
seen during banana states with the exception of 4U1728- 



4.1. The nature of 4U1915-05 

TIh p r-laQQifipatior, intmHuppH hy H^ingw aT1 H yar, r|pr Mendez and Van der Klis, 1999| ) and 4U1735-44 ( [Ford 



34 (|Ford and Van der Klis, 1998|; Strohmayer et al., 1996 



Klid (|Tq«q|) H^ngi^W twn r-la^PQ nf T.MYRq tt,P 7 et al., 1998b|). This gives further support to the idea that 



and the atoll sources. We list the main evidence sug- 
gesting that 4U1915-05 is an atoll source. First, 4U1915- 
05 has always been detected at X-ray luminosities below 
~ 1.5 x 10 37 ergs s _1 , which is typical for atoll sources. 
The Z sources are roughly 10 times brighter. Second, 
the orbital period of 4U1915-05 (50 min) is among the 
shortest observed for LMXBs. Now, it is observed that Z 
sources have longer periods (more than 10 h) than atoll 
sources, consistent with t he idea that their c ompanions 
are giants or subgiants ( Van der Klis, 1995| ). Further- 
more, almost all persistent dippers and bursters belong 
to the atoll class (Cyg X-2 and G X 17+2 are the only 
Z sources to display type I bursts ( Kahn and Grindlay. 
1984 |Smale, 1998| ; |Sztajno et al., 1986; )). Finally, no Z 
track in the color-color diagram has been observed neither 
in the GINGA data QYoshida, 1992| ), nor in the RXTE 
data. Hence, this favors the idea that 4U1915-05 belongs 
to the atoll class, as sugg ested by Yoshida ( |1992| ). 

4.2. The states of 4U1915-05 

Identifying the different atoll states in 4U1915-05 is not 
straightforward. From the hardness-intensity diagrams, 
we can see that the emission, on average, varies along the 
observations from a HS to a LS that could correspond 
to banana and island states respectively. However, we do 
not observe clearly separated regions in the color-color 
diagram but rather a continuous banana shape. Further- 
more, in its LS, 4U1915-05 does not display HFN which 



4U1915-05 was in a banana state during our observation. 
The highest intensity state where no HFQPO is detected 
(February and March observations) was therefore more 
likely the upper banana state. The lowest intensity state 
was likely the lower banana state, although no HFQPO is 
detected during May 18th to 22nd observations where the 
count rate is the lowest; maybe due to a lack of sensitivity 
(the RMS upper limits corresponding to these observa- 
tions and quoted in Table ^ are not really constraining) . 

4-3. Tracking the accretion rate with S a 

It has been shown that the relation between HFQPO fre- 
quency and intensity could be much more complex than a 
roughly one to one correlation as for example in 4U1820- 
30 (Smale et al., 1997). Numerous branches are clearly 



visible in the HFQPO frequency versus intensity diagram 
made with an extensive set of data from 4U1608-52 and 



4U1728-34 (IMcndcz et al., 1999| ; |M6ndcz and Van dci 



Klis, 1999| ). Zhang et al. (|1998| ) also showed that Aql X-l, 



observed in very distinct flux or count rate ranges, could 
display HFQPOs in the same frequency range. The ab- 
sence of a simp le correlation had p reviously been noticed 
in 4U1705-44 (|Ford et al., 1998a|) . For 4U1608-52, the 
HFQPO frequency was reported to correlate well with the 



count rate but only on timescales of hours (Mendez et al 



1999 ). Thus, the relation between HFQPO frequency and 



intensity seems to differ from source to source; sometimes 
being complex. On the contrary, the position on color- 
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color diagram yields a more universal and unique relation 
as in, e.g. 4U1608-52, 4U1728-3 4 flMendez et al, 1999| ; 
Mendez and Van dor Klis, 1999fl. Our study of 4U1915- 



05 confirms these results. The HFQPO frequency has a 
simpler relation with the position S a on the color-color 
than with the count rate (Fig. || and [l(]) . The frequency 
is well correlated with S a within its full range. 

Now, currently admitted models for HFQPOs, which 
are beat frequency models (BFM), involve a bright spot 
at a specified radius of the Keplerian disk, e.g. the sonic 
radius in the "sonic-point model" (Miller et al., 1998c) 
to produce the upper HFQPOs. As the accretion rate 
increases, the disk inner radius moves inwards until the 
innermost stable orbit is reached. Then, the HFQPO fre- 
quency is supposed to correlate with the accretion rate 
up to a saturation. Thus, the observed complex relations 
between HFQPO frequency and count rate show that the 
count rate is not as good as S a to track the accretion rate. 
This confirms the studies carried out with EXOSAT data 
which revealed the position on color-color diagrams as a 



bet ter indicator of accretion rate than intensity ( Hasinger 



andl Van der Klis, 198$ |Van der Klis et al., 1990 
et a |l-, 1990|) 



Hasinger 



However, the position on the color-color diagram may 
not be the only or the best indicator of the accretion 
rate or (and) of the timing behaviour. The properties of 
the HFQPO RMS amplitude found in 4U1915-05 (Fig. 
H and [n]) may be new evidence for this. Indeed, in the 
sonic-point model, the HFQPO amplitude is expected to 
decrease as the accretion rate increases: the optical depth 
and electron density increase near the neutron star, and 
the oscillations genera ted there are more attenuated dur- 
ing their propagation ( Miller et al., 1998c ). Therefore, one 
would expect the HFQPO RMS to be anticorrelated with 
the accretion rate. However, for 4U1915-05, the HFQPOs 
RMS does not seem to be anticorrelated with S Q but 
rather may be anticorrelated with the count rate (Fig. 
[l0| and ^). Thus, the timing properties do not probably 
depend only on the position S a . Furthermore, it is not 
clear why the HFQPOs RMS amplitude and frequency 
would be driven by different parameters. 

Moreover, we have shown that 4U1915-05 can display 
different timing behaviours (presence or absence of HFQ- 
POs) during observations overlapping in the color-color 
diagram. The same situation occurs in other sources, e.g. 
4U1636-53 QMendez, 1998| ). 

One explanation may be that these timing behaviours 
are in fact the same but seem different in appearance be- 
cause of instrumental limitations. For example, a QPO 
may be present in two observations overlapping in the 
color-color diagram but not detected in one of them be- 
cause one of the detectors has been switched off leading 
to a loss of sensitivity (see e.g. the RMS upper limits in 
Table @). 

A second explanation may be that within the times- 
pan for which S a is computed, the source might have 



moved in and out the S a span, washing out temporal 
signals. 

Another explanation is that the source has actually 
intrinsically different timing behaviours. In this case, this 
is evidence that the position in the color-color diagram is 
not the best indicator of the timing behaviour of a source. 
Then two different conclusions may be considered. First, 
assuming that the timing behaviour is governed by the 
accretion rate, then S a is not a good indicator of the 
accretion rate. This quantity may be better represented 
by a combination of parameters or by other parameters 
than S a . Second, assuming that S a is a good indicator 
of the accretion rate, then the timing behaviour is not 
governed only by the accretion rate. 



4-4- Low and high frequency quasi-periodic oscillations 

LFQPOs between a few and ~ 80 Hz have now 
been observed in several atoll sources (se e Table |ij) 



In 4U0614+09 (Ford, 1997), 4U1608-52 (Yu et al 



1997), 4U1702-43 (Markwardt et al., 199S), 4U1728-34 



( [Btrohmaycr et al., 1996 ; Ford and Van der Klis, 1998 ) 



KS1731-260 (|Wijnands and Van der Klis, 1997j ) and 
4U1735-444 flWijnands et al., 1998cj ), such LFQPOs are 
detected simultaneously with HFQPOs. In the case of 
4U0614+09, 4U1728-34 and 4U1702-43, the LFQPO fre- 
quency is observed to increase with the HFQPO fre- 
quency and with the intensity of the source. First, this 
indicates that the LFQPOs detected in atoll sources are 
comparable to the horizontal branch QPOs (HBOs) de- 
tected in Z sources, suggesting that similar physical pro- 
cesses are at work in both kind of LMXBs. Second, this 
indicates that LF and HFQPOs are produced either by 
related mechanisms or by the same mechanism seen un- 
der different aspects or occuring at different locations in 
the disk. Our study of 4U1915-05 confirms both these 
conclusions. For the first time, we show a LFQPO in an 
atoll source varying in frequency and RMS amplitude as 
a function of the position on the color-color diagram (Fig. 
|l0| ) the same way as HBOs in Z sources. We note however 
that in 4U1915-05, the LFQPOs are not strictly corre- 
lated with the source intensity, which is a characteristic 
of HBOs and also of the atoll sources mentioned above 
showing LFQPOs and HFQPOs. 4U1915-05's behaviour 
is thus different in this respect. This may be explained by 
the fact that in 4U1915-05, the intensity is not correlated 
with the position in the color-color diagram over the en- 
tire range where LFQPOs are detected. On the contrary, 
the intensity is strictly correlated with the position in the 
color-color diagram in horizontal branches ( Van der Klis J 



1995 ) and probably in some portions of the track of atoll 
sources. In 4U1915-05, we clearly show (Fig. |l0| and ||) 
that S a , and not the count rate, is the key parameter 
of the quasi-periodic variability at low frequencies. Thus, 
the HF and LFQPOs frequencies depend the same way 
upon S a - 
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Psaltis et al. (1999a) have shown that the frequency of 
the QPOs detected between 0.1 and 100 Hz in numerous 
Z, atoll or black hole binaries, simultaneously with HFQ- 
POs followed one of a small number of correlations with 
the lower HFQPO frequency. For 4U1915-05, we could de- 
tect simultaneously the LFQPO and the lower HFQPO 
only during three segments. On the diagram derived by 
Psaltis et al. ( |l999a| , see their figure 2), these points would 
fall in the region where the two main correlations merge. 
Hence, our RXTE observation does not help to conclude 
about the unified picture suggested for QPOs. 

The fact that HBOs and HFQPOs have now been ob- 



is the mean peak separation between the twin HFQPOs 
detected in 4U1915-05. Furthermore, we used the same 
parameters as Stella and Vietri ( 1998| ) for the neutron 
star: mass M — 1.97M Q , ratio I45/M — 1.98 where Z45 is 
the moment of inertia in units of 10 45 g cm 2 . With these 



served simultaneously in Z sources, as in Cyg X-2 (Wi- 



jna:ids et al., 1998a) calls into question the models pro- 
posed for QPOs. Indeed, the magnetospheric BFM was 



first proposed to explain the HBOs in Z sources (Al- 



par and Shaham, 1985). Later, the same model was sug- 
gested to explain the HFQPOs (Strohmayer et al., 1996) 
observed both in atoll and Z sources. But simultaneous 
LFQPOs and HFQPOs are difficult to explain with the 
same model. This requires that the disk is still present 
inside the radius where the magnctosphcrc cou ples the 
disk and channels the matter ( Miller et al., 1998c ). In this 
case, HBOs and HFQPOs can be expl ained b oth with the 
magnetospheric BFM. Psaltis et al. ( 1999b ) have shown 
tha t the correlations observed between HBO and HFQPO 



values, Stella and Vietri ( 1998) ) found that the precession 
frequency matched the LFQPO frequency for three dif- 
ferent atoll sources (4U1728-34, 4U0614-09 and KS1731- 
260). Since the neutron star parameters used yield preces- 
sion frequencies close to the maximum values that can be 
reached keeping the parameters in the ranges allowed in 
classical neutron star models, it is unlikely that the Lense- 
Thiring model can be pushed up to match the LFQPOs in 
4U1915-05. The precession frequencies are actually lower 
than the observed LFQPOs frequencies for any values of 
M and I45/M within their reasonable respective ranges. 
This would also be the case if the spin frequency were 
lower than the one inferred from the twin peaks separa- 
tion (348 Hz). On the opposite, the precession frequencies 
would roughly match the LFQPOs if the spin frequency 
were twice that value (~ 700 Hz). With a spin frequency 
of 348 Hz, matching the observed LFQPO frequencies 
with precession frequencies would require a ratio I45/M 
of ~ 3.5 which is too large. Indeed, this ratio is inferred 
to be in the range ~ 0.5-2 for realistic rotating neutron 



fre quencies in Z sources were indeed consistent with this 
model. However, this result has been shown so far for Z 
sources only. The striking similarities between the LFQ- 
POs detected in 4U1915-05 and HBOs suggest that both 
phenomena could be interpreted within a same model. 



star models (Stella and Vietri, 1998; Markovic and Lamb 



1998; Kalogera and Psaltis, 1999; Miller et al., 1998a) 



But! the magnetospheric BFM may not be easily extrapo- 
lated to atoli sources because of their lower interred mag- 
netic fields and accretion rates. 

Another possible interpretation comes from Stella and 
Vietri (1998) who proposed that the LFQPOs observed 



Similar conclusions have been reached for 4U1728- 
34 ( Ford and Van der Klis, 1998 , from a large r set of 
data than t he one used by Stella and Vietri ( 1998 )) 
4U1735-44 flWijnands et al., 1998c|) , 4U1702-42 flMark 



wardt et al., 1999, in this case, the required ratio I45/M 
is only 2.3 which is closer of the accepted range than for 
other sources) and for the Z sources GX 17+2, GX 5-1, 



GX 340+0, Cyg X-2 and Sco X-l (Stella and Vietri, 1998 



aro md 15-50 Hz in several atoll and Z (horizontal branch) 



Jonker et al., 1998; Psaltis et al., 1999b; Kalogera and 



sources could result from the precession of the innermost 
disk regions. In the first version of this model, the LFQPO 
frequency is the nodal precession (v no d) of slightly titled 
orbits in these regions. This precession frequency (rela- 
tivistic -Lense-Thirring- and classical contributions) de- 
pends on the equation of state of the neutron star but also 
on its spin frequency and on the Keplerian frequency of 
the innermost accretion disk region. In the framework of 
BFMs, these latter frequencies are inferred for LMXBs 
displaying twin HFQPOs, from the frequency separation 
of the peaks and the frequency of the upper peak re- 



Psaltis, 1999] ). For the atoll source 4U1608-52, the preces- 
sion frequency model does not seem to apply either: the 
spin frequency predicted from the Lense-Thirring model 
(Yu et al., 1997) is inconsistent with the one inferred from 
the frequency separation of the twin HFQPOs later ob- 
served (|Mcndcz et al., 1998b|; [Mendez et al., 1998c[). In 



summary, the LFQPOs roughly follow the dependence 
upon the Keplerian frequency expected within the Lense- 
Thirring model but occur at a frequency lower than the 
precession frequency for a growing group of binaries that 
now includes 4U1915-05. Thus, in order to explain the 



spectively. The precession frequency is predicted to vary esis are required ( 


Markovic and Lamb, 1998|; 


Morsink and 


app 
free 


roximately as the square of the innermost Keplerian Stella, 199E ; 


Armitage and Natarajan, 1999; 


Kalogera] 


uency. For 4U1915-05, the LFQPO frequency varies as and Psaltis, 1999; Schaab and Weigcl, 1999|). 


Another 



Vjf 1 where vk is the inferred Keplerian frequency (see Fig. 
|l3| ). This is not far from the expected quadratic depen- 
dence within the errors. However, the observed LFQPO 
frequency is ~ 2 times greater than the expected pre- 
cession frequency drawn as a dashed line on Fig. [l3| We 
assumed a neutron star spin frequency of 348 Hz which 



possibility is that the LFQPOs correspond to the second 
harmonics of the precession frequency 2v no( i instead of 
v n od as proposed in a second version of the Lense-Thirring 
model, maybe because of a modulation at twice the pre- 
cession frequency generated at the two points where the 
inclined orbit of the blobs intersects the disk ( Vietri and 
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Stella 1998; [Morsink and Stella, 
199|; |Stella et al., 1999| ). 



1999: Stella and Vietri 



However, Psaltis et al. ( |l999b[ ) and Schaab and Weigel 



(1999) noticed another discrepancy with the model: the 
LFQPO frequency does not follow the expected quadratic 
dependence on the Keplerian frequency when the latter 
frequency is greater than 850 Hz in Z sources. In addi- 
tion, radiation forces in the disk could significantly affect 
the precession frequencies and thus possibly challenge the 
Lense-Thirrin g model as the explanation for the LFQPOs 
flMillcr, 19991) . Thus, the LFQPOs origins and relations 



with HFQPOs remain unclear. 

Another model proposes to interpret both the LF and 
HF QPOs within the same framework: the two-oscillator 



model (Oshcrovich and Titarchuk, 1999b 




Titarchuk 


and 


Oshcrovich, 1999]: 


Titarchuk et al., 1999 


). In this 



model, the lower kHz QPO is the Keplerian frequency 
(vk ~ 700 Hz) at the outer edge of the boundary layer 
between the Keplerian disk and the neutron star. Blobs 
thrown out from this region into the magetosphere os- 
cillate in both radial and perpendicular modes to pro- 
duce an upper kHz QPO (at Vh ~ 1000 Hz ) and a low 
frequency QPO (at vl ~ 50 Hz) respectively. Further- 
more, viscous and diffusive processes within the bound- 
ary layer produce respectively a low frequency QPO (at 
vy ~ 30 Hz) and a break (at v\, ~ 8 Hz) in the PDS. In 
this model, the angle S between the rotational frequency 
v = Vt/2ir of the magnetosphere and the normal to the 
neutron star disk is not zero. 

The various observed QPOs of Sco X-l, 4U1728-34 
and 4U1702-42 have been succesfully identified in the 



framework of the two-oscillator model ( Titarchuk et al.,| 
199|; psherovich and Titarchuk, 1999a| ). The angle 8 de- 
rived is 8.3° ±1.0°, 5.5° ±0.5° and 3.9° ±0.2° respectively 
for each source. Furthermore, the predicted variation of 
Vb as Vy 6 has been checked for number of atoll and Z 
sources ( |Titarchuk et al., 1999| ). 

For 4U1915-05, the identification of the observed 
QPOs is not straightforward. 4U1915-05 does not show 
any break in its PDS. Only one LFQPO is detected in the 
frequency range where QPOs are predicted to be present 
at vy, vl and 2vl- 

We tried however to test the two-oscillator model, ft, 

(fl/ir) 2 (see 



vk and Vh follow the relation v\ 



e.g., equation 2 in Osherovich and Titarchuk (1999a)), so 
that the rotational frequency can be derived when twin 
HFQPOs are detected simultaneously and interpreted as 
vk and Vh- 

The bottom panel of Fig. |lj shows the inferred ro- 
tational frequency as a function of the lower HFQPO 
frequency, for the four segments where twin HFQPOs 
are detected above 550 Hz (see Table 0). Theoretically, 
ft depends on the magnetic structure of the neutron 
star's magnetosphere. In the case of 4U1915-05, there 
are not enough simultaneous detections of vk and Vh to 
reconstruct the ft profile. We thus use the approxima- 
tion ft = Oq = const, as was done also for 4U1702-42 



zl> 
450 - 



300^ 
500 



600 650 700 

Lower HFQPO frequency (Hz) 



Fig. 14. Inferred rotational frequency (bottom panel) 
and angle 6 (top panel) as a function of the lower HFQPO 
frequency within the two-oscillator model. The dashed 
lines represent the mean values. The four points in the 
bottom panel are inferred from the 4 detections of twin 
HFQPOs above 550 Hz. The three points in the top panel 
correspond to the three cases where a LFQPO is simul- 
taneously detected. 



( |Osherovich and Titarchuk, 1999a ). As shown in Fig. |lj 
(bottom panel), the rotational frequency is indeed consis- 
tent with being constant with v = Qq/2tt = 373 ± 36 Hz. 
We note that the pair of QPOs at 224 and 514 Hz would 
correspond to a rotational frequency of 231.2 ± 1.3 Hz 
inconsistent with the previous value and not considered 
here. 

Knowing ft, the angle 6 can be derived when vk and 
vl are measured simultaneously (see e.g. equation 5 in 
Osherovich and Titarchuk (1999a)). The top panel of Fig. 
|l4| shows the inferred angle 5 for 4U1915-05, assuming 
that the LFQPO frequency is vl- S is consistent with 
being constant with a mean value S = 4.6° ± 3.1°. It is in 
the range of values found for the other sources, but the 
error is particularly large. 

We note however that the relation between ft, vk and 
Vh written above requires v% > 2v, which becomes Vh > 
746 Hz in the case of 4U1915-05. Now, according to Fig. 
|Io| (left) showing the HFQPO frequency as a function 
of S a , when only one HFQPO is detected, it seems to 
be the upper HFQPO. If we assume that it is actually 
the case and that its frequency is Vh, then the observed 
frequency range of v^ would begin at ~ 500 Hz. This 
would be inconsistent with the condition v^ > 746 Hz 
derived according to the assumptions considered above. 
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4-5. Implications for the neutron star in 4V 1915-05 

Observations of HFQPOs allow to derive constraints on 
the neutron star present in the system. In the framework 
of beat frequency models, the frequency difference be- 
tween twin HFQPOs is interpreted as the spin frequency 
of the neutron star. For 4U1915-05, the mean frequency 
separation is 348 Hz for four segments of observations. 
This is consistent with values reported in other LMXBs 
and would imply a neutron star rotating with a 2.8 ms 
period. However, a frequency separation of 290 Hz is de- 
tected for a pair of QPOs at 514 and 224 Hz. The spin 
frequency detection at 348 Hz needs to be firmly con- 
firmed. 

5. Conclusions 

We have analyzed the 1996 RXTE data of the X-ray 
burster and dipper 4U1915-05. We confirm that it is an 
atoll source that was probably in its lower and upper ba- 
nana branch during the RXTE observations. 

We report the discovery of both LFQPOs (5-80 Hz) 
and HFQPOs (200-1300 Hz) and weak evidence for a 
0.2 Hz QPO. We show a correlation between the LFQ- 
POs and the HFQPOs frequencies suggesting that both 
QPOs are produced either by related mechanisms or by 
the same mechanism seen under different aspects or oc- 
curing at different locations in the disk. Furthermore, 
the LF and HFQPOs frequencies both positively corre- 
late with the position S a of the source in the color-color 
diagram that is a key parameter of the timing behaviour 
though probably not the only one. 

Four twin HFQPOs are detected above 550 Hz with 
a frequency separation consistent with being constant 
(348 ±12 Hz) suggesting a spin frequency of 2.8 ms for 
the neutron star. A pair of QPOs below 550 Hz is also 
detected with a frequency separation of 290 Hz. 

Many more observations have been performed by 
RXTE. They should help to confirm the picture of the 
aperiodic variability discovered from 4U1915-05. Espe- 
cially, the 0.2 Hz QPO and signatures of the neutron star 
rotation should be searched for in the persistent, dipping 
and bursting emissions. 
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